arXiv: 1506.01471 vl [math.NA] 4 Jun 2015 


Minimal Gersgorin tensor eigenvalue inclusion set and its numerical 

approximation 


Chaoqian Li and Yaotang Li 

School of Mathematics and Statistics, Yunnan University, Kunming, Yunnan, P. R. China 650091 
Email: lichaoqian@ynu.edu.cn, liyaotang@ynu.edu.cn 

Abstract 

For a complex tensor A, Minimal Gersgorin tensor eigenvalue inclusion set of A is 
presented, and its sufficient and necessary condition is given. Furthermore, we study its 
boundary by the spectrums of the equimodular set and the extended equimodular set for 
A. Lastly, for an irreducible tensor, a numerical approximation to Minimal Gersgorin 
tensor eigenvalue inclusion set is given. 
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1. Introduction 

Let N = {1,2,... ,n}. For a complex (real) order m dimension n tensor A = ( 0 ^ 1 ...^^) 
(written A G (]R["*’"1), respectively), where 

G C ij = 1,. ■. ,n, j = 1,... ,m. 

we call a complex number A an eigenvalue of A and a nonzero complex vector x an eigenvector 
of A associated with A, if 

= Axt"*-!!, (1.1) 

where Ax"^~^ and are vectors, whose zth components are 

(Ax ')i — ^ ^ tlii2...i^Xi2 * * * 

and 

(x^^-% = xT-\ 

respectively. Note that there are other definitions of eigenvalue and eigenvectors, such as, 
H-eigenvalue, D-eigenvalue and Z-eigenvalue; see [!3 M [19]. Obviously, the definition of 
eigenvalue for matrices follows from the case m = 2. 

Tensor eigenvalues and eigenvectors have received much attention recently in the literatures 
[iiiaiiiiToKiaiinKiaiiTiiia. Many important results on the eigenvalue problem of matrices 
have been successfully extended to higher order tensors; see [miiiiiKiniiiniiiiiMiiii]. 
In |15] , Qi generalized Gersgorin eigenvalue inclusion theorem from matrices to real supersym¬ 
metric tensors, which can be easily extended to generic tensors; see m- 
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Theorem 1.1. [15] Let A = and a(A) be the spectrum of A, that is, 

a{A) = {A e C : Ax'^-^ = x £ C^VIO}}. 


Then 



( 1 . 2 ) 


ieN 


where 



and 



It is easy to see that when m = 2, Theorem 11.11 reduces to the well known Gersgorin 


eigenvalue inclusion theorem of matrices 13111]. Here, we call Ti (^) the i-th Gersgorin tensor 


eigenvalue inclusion set. Note that ri(Al) is a closed set in the complex plane C, Hence, r(./l), 
which consists of the n sets ri(.4,), is also closed and bounded in C. 

Also in Qi obtain another interesting result on Gersgorin tensor eigenvalue inclusion 
set r(A). 

Theorem 1.2. [T5] IfTi(A) is disjoint with the other rj(Al), j ^ i, then there are exactly 
(m — 1)”“^ eigenvalues which lie in ri(A). Furthermore, if all of ri(A), i = li,l 2 , ■ ■ ■, Ik o,re 
connected but disjoint with the otherVj^A), j ^ i, then there are exactly fc(m—1)"“^ eigenvalues 
which lie in U r,(A). 


i—h 


In this paper, we also focus on Gersgorin tensor eigenvalue inclusion set. And we present 
Minimal Gersgorin tensor eigenvalue inclusion set, give its sufficient and necessary condition, 
and research its boundary. For an irreducible tensor, we give a set which approximates to 
Minimal Gersgorin tensor eigenvalue inclusion set. 


2. Minimal Gersgorin tensor eigenvalue inclusion set 


In this section, we present Minimal Gersgorin tensor eigenvalue inclusion set and study its 
characteristic. First, a lemma is given. 

Lemma 2.1. [53] Let A = and D = diag{di,d 2 , ■ ■ ■ ,dn) be a diagonal 

nonsingular matrix. If 


m — 1 


B = {b,,...,J = DD- -D 


where 



then A, B have the same eigenvalues. 


From Lemma im we obtain the following tensor eigenvalue inclusion set. 
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Theorem 2.1. Let A = and x = (xi,X 2 , ■ ■ ■ be an entrywise positive 

vector, i.e., x = {xi,X 2 , • ■ ■, XnY' > 0- Then 

a(. 4 )cr^(^)= U rr(. 4 ), ( 2 . 1 ) 

ieN 


where 

rr(^) = 

Furthermore, 


zeC: |z-a....,| <r“(.4)= ^ 




xT~^ 


r(^) c f| 


( 2 . 2 ) 


£C>0 


m—1 

Proof. Let X = diag{xi,X2, ■ ■ ■, Xn) and B = XX ■ ■ ■ X = It is obvious 

that X is nonsingular and cr{B) = criA) from Lemma l2.II From Theorem ll.il we have 

a(S)cr(s)= Ur*( 6 ), 

ieN 


where 

r,{B)=lz€C:\z -< r,(6) = ^ | . 

i2. im€N, 

Note that bi...i = at—i and ri{B) = rf{A). Hence, T'i{B) = rC*(.4), consequently, r(H) = 
r’'*(.4). Therefore, 

a{A) =cr{B) cr{B) = r^^{A). 

Furthermore, for any x > 0, x G (EH) also holds. Hence, (1^ follows. 

The set p] F’’* {A) in Theorem 12.11 is of interest theoretically because it provides a set, 

3;>0 

containing all the eigenvalues of A, which is called Minimal Gersgorin tensor eigenvalue inclusion 
set defined as follows. 

Definition 2.1. Let A = (uij...^^) G Cf"*’"!. Then 

r^(A)= f| r^^(A) (2.3) 

x>0 

is Minimal Gersgorin tensor eigenvalue inclusion set of A, related to the collection of all 
weighted sums, rf(A), where x = {xi,X 2 , ■ ■ ■, x„)^ > 0. 

Next, a sufficient and necessary condition for the elements belonging to Minimal Gersgorin 
tensor eigenvalue inclusion set is provided. Before that, we give some results involving the 
Perron — Frobenius theory of nonnegative tensors |23] . Given a real tensor A, we call A 
nonnegative, denoted by M > 0, if every of its entries is nonnegative. 
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Definition 2.2. |IJ Definition 2.1] A tensor A = S is called reducible, if there 

exists a nonempty proper index subset I C N such that 

l^zi ■ ■ ■ , 2^71 1 — ^ h, \fi2, . • . , im ^ h. 

If A is not reducible, then we call A irreducible. 

Lemma 2.2. [53] If A G is nonnegative, then the spectral radius 

p{A) — max{|A| : A S cr(Al)} 

is an eigenvalue with a nonnegative eigenvector x 0 corresponding to it. Moreover, if A is 
irreducible, then p{A) > 0 and x is positive. 

Lemma 2.3. |23l Theorem 5.3] Let A G be nonnegative. Then 


p{A) = max min 




x>0,x^0xi>0 ^ 


= mm max ■ 

x>0,x^0xi>0 x^- 


{Ax^-^fi 


Remark 2.1. The first equality is only given in Theorem 5.3 of IMi- And similar to the proof 
of Theorem 5.3 of ^23] . the second is proved easily. 


From Lemma l2.31 we can obtain the following result. 
Corollary 2.1. Let A & RI"*’”! be nonnegative. Then 


p{A) = sup rriin 






= inf max ■ 


m—1 


c>0 X^ a;>0 ieN x.^ 

From Corollary 12.11 we have the following result. 

Lemma 2.4. Let A = G and z be any complex number. And let B{z) = 

G R[™’"1, where 


bi...i — jz 


= a*' 


for i G N and Si, 


= 0 . 


Then B{z) possesses a real eigenvalue v(z) which has the property that if X G a{B{z)), then 

Re{X) < v{z). 


Furthermore, 

(2.4) 

^ ' x>o ieN 

Proof. Let p = max \z — ai...i\, and C = (c^^... G RI"*’”!, where 

i^N 

— IJ> \z di-.-i I ClTld, - -- ,im. I J ^ ^7 ^ii2 ■ ■ ■ (^■^) 

Then B{z) = —pI + C, where C > 0. From Equality (II.IL it is obvious that A(C) G a{C) 
if and only if —p + A(C) G a{B{z)). Moreover, From Lemma 15^ p{C) is an eigenvalue of C 
with a nonnegative eigenvector. Hence, —p + p{C) is an eigenvalue of B{z) with a nonnegative 
eigenvector. Let 


v{z) = -p + p{C). 


(2.6) 
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Obviously, v(z) G a{B{z)) is real. And for any A e a(B(z)), we have 

Re{X) = Re{X + ^ — fi) = Re{X + /i) — /r < |A + /r| — /i < p{C) — p = v{z). 
Furthermore, by Corollary 12.11 we have 


p(C) = inf max ■ , 


Then 


v{z) = -p + piC) 

= —p + mt max- 


= inf max 
a;>0 ieN 


= inf max 


x>0 ieN 

(i-pX + C)x^-^\ 




a;>0 i^N x 


771—1 


The conclusion follows. 

For v{z) in Lemma l2.41 we give the following property. 

Proposition 2.1. Let A = G and v{z) be defined as Lemma Then v{z) is 

uniformly continuous on C. 

Proof. Let B{z) G rI™’”! be defined as Lemma 12.41 Then, similar to the proof of Theorem 

[121 


v(z) = inf max 
£c>o ieN 


iB{z)x'^-^fi 


iiif max{rf (A) - \z 

x>0 i^N 


Note that for any z, z gC and a; > 0, x G R", 




maxjrf (A) - \z - ai...^\} = 

ieN 

> 


maxfrf (A) - \z - z + z - at...A} 

iGN 

max{rf (A) - i\z - 5| + |z - a^...i|)} 
max{rf (A) - |z - ai...*|} - |z - z|, 

iGN 


that is, 


Hence, 


max{rf{A) — \z — ai...i\} < \z — z\ + inax{rf{A) 


\z - 


inf niax{rf (A) - |z - a^...* |} < 

x>0 t^N 


inf Uz-S\+ max{rf (A) - \z - 

a:>0 t^N 

\z-z\+ inf inax{rf (A) - |z - a*...*|}, 

x>o leN 


which implies 
Similarly, we can obtain 


v{z) — v{z) < |z — z|. 
v{z) — v(z) < |z — z|. 
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Therefore, 

\v{z) — v[z)\ < |z — 5|. 

This implies that v{z) is uniformly continuous on C. 

We now establish the sufficient and necessary condition for r^(^). 

Theorem 2.2. Let A, B{z) and v{z) be defined as Lemma lK^ Then 

z € r'^(.4) if and only if v{z) > 0. (2-7) 

Proof. Assume that z G r'^(A). From Definition 12.11 we have that for each vector a; > 0, 
z G T’' (A). Hence, there is an io G N such that 


equivalently, 

Note that for any i G N, 


k-a*o-iol < rl(A), 
rfoi-A.) - \z - \ >0. 


Then - —— > 0, which implies that for each vector a:: > 0, 

^*0 

(S(z)a:—1), ^ „ 
max-- > 0. 


By Lemma [2.41 


ieN X. 


'j{z) = inf max 


,m— 1 


(S(z)a™-i), 


> 0 . 


x>o ieN x" 

Conversely, suppose that v{z) > 0. From Equality (12.4L then for each vector x > 0, there 
is an i G such that 

0 < v{z) < = r^A) -\z- a....,|. 

Hence, \z — ai...i\ < rf{A), and then z G Vf {A) which implies z G r’'''(A). Since this inclusion 
holds for each vector x > 0, we have z G F'^(A) from Definition 12.11 


3. Boundary of r^(Al) 

In Section 2, a sufficient and necessary condition for the elements belonging to Minimal 
Gersgorin tensor eigenvalue inclusion set F^(M), is given. We in this section focus on the 
boundary of F^(M), and establish relationships between its boundary (see Definition Id.Ill , the 
spectrum of the equimodular set for A (see Definition 13.21) . the spectrum of the extended 
equimodular set for A (see Definition 13.31) and r^(M). 

Definition 3.1. Let A G The boundary of Minimal Gersgorin tensor eigenvalue inclu¬ 

sion set F^(A), denoted by 9r^(M), is defined by 

dT^{A) = T%A) Pi (r«(A))', 

where (r'^(A))' is the complement ofT^{A), and F^(A) is the closure o/F^(M). 
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Definition 3.2. Let A = The equimodular set of A, denoted by is 

defined as 

~ {_Q ~ {Qil-'-im) ^ ^ ■ Qi---i ~ 1^222 ~ l^2^2---*m 17 i ^ LI ^ ^2^2 ■■■^m “ 

Definition 3.3. Let A = (aij...i^) € Cl™’"!. The extended equimodular set of A, denoted by 
r2(.4,), is defined as 

“ {_Q ~ ^ ^ ~ l^222-"2Tn I — I ^222-"2m I: ^ ^ ^222" "2m “ 

A sufficient and necessary condition for the point lying on the boundary of r^(A) is given 
as follows. 

Proposition 3.1. Let A and v{z) be defined as Lemma \K^ z € cir^(A) if and only if v{z) = 
0, and there is a sequence of {zj}fLi C (r^(A))' (i.e., v(zj) < 0 for all j > 1) such that 
lira Zj = z. 

j^oo 

Proof. Note that r^(A) is a compact set in the complex plane. Hence, from Definition 13.11 

ar«(A) = W(A) Pi (r«(A))' = r^(A) p| (r«(A))'. 

Furthermore, by Theorem 12.21 we have that 

z G (r^(A))' if and only if v{z) < 0. (3.1) 

Therefore, if z G 9r^(A), then, from (12.7|) and (13.IL v{z) = 0. Note that (r^(A))' is open and 
unbounded, then, by (13.11) and z G (r^(A))', there is a sequence of {zj}jTi C (r'^(.4))' such 
that lim Zj = z, where v(zj) < 0 for all j > 1. 

j-yoo 

Conversely, it is obvious by assumption that z G (r^(A))' and z G r'^(A) = r^(A), that is, 

z G 5r^(A). 


The proof is completed. 

Proposition 3.2. For any A = (ai^...i„) G and any z G C with v{z) = 0, there is a 

tensor Q = {qi-^-.-i^) G D(A) for which z is an eigenvalue of Q. Then 

dT^iA) C a{n{A)) C a(D(A)) C T^{A), 

where a{n{A)) = IJ cr(I?) and a {Cl (A)) = IJ cr(X>). 

'DGn{A) -DeniA) 

Proof. From Definition 12.11 [3?^ and [3Al it is obvious that D(A) C Cl{A) and cr(D(A)) C 
ct(D(A)) C r^(A). Hence, we next only prove 9r^(A) C tT(D(A)). 

If z G C is such that w(z) = 0, then, from the proof of Lemma 12.41 there is a vector 
y = (j/i,... ,yn)'^ > 0, y 7 ^ 0 such that B{z)y = 0, where S(z) is defined as Lemma This 
implies that for any k € N 

'y ^ I®fcfc2"'fcm 12/^2 2/fcm • 

k2 . km€N, 


\z-ak...k\y^ ^ 


(3.2) 









Now, let ilik satisfy 


z - ak-k = \z- afc...fc|e*'^^ fc G iV, 
and Q gC[™'" 1, where 

Qk‘-k — ^k---k Oyfld Qkk2‘“km — \^kk2---km.\^'^^ ^ ^ ^ ^kk2---km — 0 . 

Hence, from Definition 13.21 Q G n(^). Moreover, By considering the k-th entry {Qy"^~^)k of 
Qyin-i^ we have that for any k € N, 


{Qy )k — ^ ^ Qkk2---kmyk2 ' ' ' TJkm 

k2,...,krr,eN 

/ 

= Qk-kyk + e 'y ^ \^kk2-"km\yk2 ' ' ' yk„ 


— (z — (z — ak---k))y^ + 


m-l I ii/’fe 




fc 2 ,...,fcmeAf, 
Y ^kk2---kjn=^ 


\ 




/ 

V 


= ^J/r ^ + -|^ - afc...fc|yr Wkk2---kjyk2---yk„ 


k2,---,km€N, 

^kk2 ■■■km 


= zj/™ \ {By Equality (1??^ ) 


that is, 

Note that y ^ 0, then z is an eigenvalue of Q G which shows that v{z) = 0 implies 

z G (T(n(^)). From Proposition 13.11 we have that for each point z G 9r^(^), v{z) = 0, 
consequently, z G cr(n(^)). Hence, 9r'^(^) C a{V,{A)). The proof is completed. 

For the sets cr(H(^)) and F^(yl), we have the following result. 

Proposition 3.3. Let .4 = ( 0 ^^...^^) G Ct™’"!. Then 

ain{A)) = r«(M). 

Proof. From Proposition 13.21 we have (t(H(. 4)) C r^(M). Hence, we only prove F^(M) C 

a{n{A)). 

Let z G r^(.4). Then, from Theorem 12.21 v{z) > 0. And from the proof of Lemma [2.41 
there is a vector y = {yi,..., y„)^ > 0, y 0 such that 

S(z)y™-i = u(z)y[™-h, 

where B{z) is defined as LemmaHence for any k G N, 

(^\z Q/e.../e| “h ^ ^ \^kk2---km 12/^2 ’ ' ’ Vkyn 

^2 ’ •' • >^Tn ^ ^) 

^kk2--km 


(3.3) 
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Now, let Q = G with 

Qk---k — ^k---k OjTid Qkk2---krn — f-^k^kk2---km : ^ ^ -^5 ^kk2---km — 0 , 

where 


/ 


y . \^kk2 • ■ kjyi \yk2 "'Vkj, 

^kk2 -km 


\ 


-«(z)yr ^ 


A^fc = < 


^ ^ I fc2 • • • AItti I 2/fc2 ' ' ' 2/^77 

k2,---,km€N, 

^kkn-km~^ 


1 , 


*/ E \akk2 -km\yk2-■-ykm > 

k2,...,km^N, 

^kk2 - ■ ■ krn ^ 

*/ E \a-kk2-km\yk2-■-yk^ =0. 

k2,...,km^N, 

^kkn-km.—^ 


Furthermore, from Equality (13.31) and the fact that both \z — ah...i.\v^ ^ > 0 and v{z)y]^ ^>0 
hold for any k G N, we easily obtain 0 < < 1 for any k G N. Hence, 

Q G niA). 

For the tensor Q, we have from Equality (13.311 that for any k G N, 

\z - qk-k\yT~^ -- 


z ak---k\y}^ 

( 


\ 


E \0'kk2---k^\yk2 ■ •-yk^ 


k2,...,km€^, 




= y-k 


E \0'kk2-km\yk2- ■ -ykr, 

k2,...,kmeN, 

\ <5fefe2- fcm 

^ ^ \Qkk2---km.\y!i2 ‘ ' yf^m^ 


- viz)y^ 

\ 

/ 


m—1 


k2,...,km€N, 


i.e.. 


k - gfc-fe|j/r ^ = E \^kk 2 -km\yk 2 ■■■ykrn for any k G N. 


^2 ’ •' • ykm > 

^kk2--km 


Now, the above expression is exactly of the form of Equality (liO) in the proof of Proposition 
O Hence, similar to the proof of Proposition 13.21 we have that there is a tensor V G ft{Q) 
such that 0 G cr('P). Note that V G H(Q) and Q G H(.4). Therefore, V G H(^), consequently, 
z G cr(H(^)) and r^(^) C tT(H(^)). The proof is completed. 

From Propositions 13.21 and l!01 we can obtain the following relationships. 

Theorem 3.1. Let A G Cl™’"!. Then 

dr^{A) C a{n{A)) C r^(.4). 

Remark 3.1. From Provosition \3.‘A we known that CT(r2(.4)) completely fills out F^(^), that 
is, cr(H(.4.)) = r^(.A). And from Theorem roi we known that if a {LI (A)) is a proper subset of 
F^(.4) with tT(H(^)) yA 9r^(yl), then cr(r2(yl)) necessarily have internal boundaries in F^(./l). 
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4. A numerical approximation of Minimal Gersgorin tensor 

eigenvalue inclusion set 


Unlike Gersgorin tensor eigenvalue inclusion set r(^), or F’''” {A) of (I2.1L Minimal Gersgorin 
tensor eigenvalue inclusion set r^(^) of a complex tensor A is not easy to determine numerically 
generally. In this section, for an irreducible tensor A, we give a numerical approximation of 
r^(4.), which contains r^(^). We now give a lemma, which will be used below. 

Lemma 4.1. Let A = be irreducible, and v{z) be defined as Lemma 

Then for each i G N, 

v{ai...i) > 0 . 

Furthermore, for each ai...i and for each real 6 with 0 < 6 < 27r, let > 0 be the smallest 
7 > 0 for which 


v{a,...i + 7 i( 6 ')e*®) = 0 , 
and there is a sequence of with 

lim Zj = Oi-.-i + 7 i( 0 )e*® such that v{zj) < 0 for all j > 1 . 


(4.1) 


Then, the complex interval \ai...i + te'^], for 0 < t < 7 i( 0 ), is contained in r^(4.) for each real 
0 € [0,27r), that is, the set 

(J [oi-.-i + i & N, 

e e[o,2ii-) 

is a subset ofT^{A). 

Proof. Let B{z) = G RI™’”! be defined as Lemma where. 


bi-.-i — \z and bn.^-— ,im — 

And let C = (cii... e RI'"’’"!, where 

Ci-.-i — fa \z ct^...j| and ca^-.. 7 ^ — |u, 


■ 7m 17 i G N, Si'. 


= 0 . 


|, i G Lh, Sii2... 7 ^ — 0, 


pL = maxlz — ai-.-A. Then Biz) = —pi + C. Since A is irreducible. Biz), also C, is irreducible. 

i^N 

Hence, C is an irreducible nonnegative tensor. From Lemma l2.2l there is a positive eigenvector 
of C corresponding to p{C). Therefore, similar to the proof of l2.41 there is a positive eigenvector 
of B{z) corresponding to v{z). Moreover, by Equality (12.411 . that is, 


}{z) = inf max 


{B{z)x^-^), 

^m — 1 ’ 


x>0 i£N X 

we get that for any z, there exists a positive vector y such that for all j G N, 


v(z) = 


(B(z)y^-^), 


y?-" 


(4.2) 


Then for any i G N, take z = 0 ^...^ and let x = [xi, ..., Xn)"^ > 0 be such that for any j G N , 
v{ai...i) = )j ^ particular, 

'^3 

{B{z)x'^~^)i xr A\ \ I xr A\ 

v{ai...i) — m-i ~ U (■^) Wi-'-i ai...i\ — Tj (A). 
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Since A is irreducible, we have rf(A) > 0, consequently, v{ai...i) > 0. 

Furthermore, for each 9 with 0 < 0 < 27r, + te^, for all t > 0, is the semi-infinite 

complex line. Obviously, the function v(ai,,,i + te^) is positive at t = 0, is continuous on this 
line, and is negative outside r'^(^) from Theorem 12.21 Hence, there is a smallest ^i{9) > 0 
satisfying (14.11) . And by Proposition 13.11 we get that ai,,,i + € 9r^(yl). This implies 

that for each real 0, + fy(0)e®®] is a subset of r^(A). The conclusion follows. 

We now give the following procedure for approximating Minimal Gersgorin tensor eigenvalue 
inclusion set r^(M) for an irreducible tensor A = S 

Procedure of numerical approximation 
Step 1. determine the positive numbers 

Step 2. determine the largest 7j(0) for 0 < 0 < 27r, such that aj...j -\- ^j{9)e^ € 5r^(M), 

i.e., 

v{aj...j + 7j(0)e*®) = 0, with v{aj...j + (7j(0) + £)e*®) < 0 (4.3) 

for all sufficiently small e > 0; 

Step 3. take the m points = aj...j -b 7j(0)e®® G 5r^(M) for kj^ = l,2,...,m, and 

m 

detemine the set n r (A) approximating to r^(M), where 

^30=1 

{A) = [J {z G C : \z - ai...i\ < \wkj^ - 

i^N 

Remark 4.1. (i) To determine v{aj...j), we use Equality with z = i.e., 


v{aj-j) = -^l + p{C), 

where p = m.ax.\aj...j — ai...i\ and p(C) is an eigenvalue of the irreducible nonnegative tensor 

iGN 

C defined as i2.5\) (note that the irreducibility of C is deduced by that of A), and the following 
method for calculating p{C) (see \1A \W(). i.e., if for an vector >0, let'D = C + hi, where 
h> Q, and let y^^'> = V ^, 


n(l) 




, , , = V 

(2) ^ ^ y(2) ^ p (^(2)) 


m—1 


(fc+1) 


^ V 1 > - 




and let 


Aj. = min 


(j/W). _ (y«). 

1 ’ / /I \\ m—1 ’ 5 5 5 


(ft') 






Ai — A 2 — * * * — p(A^') — p(^) H“ ^ A ''' A -^2 A 


then 
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moreover, 

lim Aj, = p(I?) =/ 9 (C) += lim Afe. (4-4) 

k —^■oo k—>00 

Hence, we can obtain convergent upper and lower estimates of v(aj...j), which do not need great 
accuracy for graphing purpose, as Example EH] shows. 

(ii) The numerical estimation of From Lemma |^.1[ there is ■jj{0) such that i4-3\ ) 

holds. Now, let z = aj...j and z = aj...j + 7 j( 0 )e*®, we have from Provosition \2.1\ that 

lj{0) > v{aj...j) > 0. 

Hence, v(aj...j + v{aj...j)e’^) > 0. Ifv{aj...j + v(aj...j)e^) = 0, then take 

7 j( 6 ») = v{aj...j) 

for which v{aj...j) can be determined by the method of (i), otherwise, v{aj...j + v{aj...j)e’'^) > 0, 
then we increase the number v{aj...j) to v{aj...j)+N, A > 0, until v{aj...j + {v{aj...j)+A)e’'^) < 0, 
and apply a bisection search to the interval [v{aj...j), w(aj...j) + A] to determine jj{0) satisfying 
O- Note here that estimates of ■jj (9) also do not need great accuracy for graphing purpose. 

(Hi) It is obvious that Wk,j^ in Step 3, kjg = 1,2,... ,m, are not only boundary points of 
r^(7l), but boundary points ofT [A), and that 

r^(^) C (7l) 

which shows that the larger m is, the better ( 7 I) approximates to r^(7l). 

Example 4.1. Consider the irreducible tensor 

.4= [A(l,:,:),Al(2,:,:),Al(3,:,:)] GCP.31, 

where 

/ 200 \ / 000 \ / 110 \ 

A(l,:,:)= 0 0 1 ,A(2,:,:)= 0 2 0 , A(3,= 0 10. 

Vool/ Vloo/ Vool/ 

We next give a numerical approximation to r^(7l). By the part (i) of Remark \4.1\ we compute 
viam) for i = 1,2, 3, and get 

?;(aiii) = v{a222) = 1.620 1 98 03, ^( 0333 ) = 1.43720383. 

Furthermore, based on the entries am = 0222 = 2 and 0333 = 1, we look for six points 

wi = am+7i(0), W2 = a333 - 73(7r), 

W3 = am + * 7i(f), W4 = am - 2 7i(^), 

W 5 = 0333 + i 73 (f), m = 0333 “ * 

of dr^{A), which are found by the method proposed in the part (ii) of Remark \4.1[ that is, 

wi = 3.62019802, W2 = -0.43720383, W3 = 2 + 21.86790935, 

W4 = 2 - 21.86790935, ws = 1 + 21.81661895, we = 1 - 21.81661895. 

Andnowr(A), r“i(^), r“=(.4), r’"^(.4) (r’"^(.4) = r’""(.4)) andr-^^iA) (r’"^(.4) =r'"«(^)) 

are given by Figures. 1, 2, 3, 4 onrf 5, respectively. In Figure 6, the set n r™'=( 7 i) , 

\fc=l / 
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which approximates to is shown with the inner boundary, and the boundary o/r(^) 

is shown with the outside. The six points are plotted with asterisks. As we can see, 

( Pi J C r(^), that is, the set which approximates to Minimal Gersgorin tensor eigen- 

k=i J 

value inclusion set is also contained in Gersgorin tensor eigenvalue inclusion set. Also, it is 
easy to see that more points of dr^{A) are given, the better {A) approximates to r^(^). 

5. Conclusions 

In this paper, we present Minimal Gersgorin tensor eigenvalue inclusion set r^(^), give 
a sufficient and necessary condition for r^(^) by using the Perron — Frobenius theory of 
nonnegative tensors, and establish the relationships between 9r^(M), (T(r2(^)), cr(r2(^)) and 
T^iA), i.e., 

5r«(^) C ain{A)) C a(f7(M)) = r^(^). 

These results obtained are generalizations of the corresponding results of matrices m to higher 
order tensors. In m, Li et al. provided two new eigenvalue inclusion sets which are contained 
in Gersgorin eigenvalue inclusion set for tensors. An interesting problem arises: what’s the 
relationship between Minimal Gersgorin tensor eigenvalue inclusion set and the sets in m ? In 
the future, we will research this problem. 
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